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43 I ABSTRACT 

HH 223 is the optical counterpart of a larger scale H2 outflow, driven by the protostellar 
Q I source VLA 2A, in L723. Its poorly collimated and rather chaotic morphology suggested the 

3-h . Integral Field Spectroscopy (IFS) as an appropriate option to map the emission for deriving 

the physical conditions and the kinematics. Here we present new results based on the IFS 
observations made with the INTEGRAL system at the WHT. The brightest knots of HH 223 
(~ 16 arcsec, ~ 0.02 pc at a distance of 300 pc) were mapped with a single pointing in the 
■ spectral range 6200-7700 A. We obtained the emission-line intensity maps for Ha, [N II] 6584 

^ \ A and [S II] 6716, 6731 A, and explored the distribution of the excitation and electron density 

Q\ ■ from [Nll]/Ha, [Sll]/Ha, and [S II] 6716/6731 line-ratio maps. Maps of the radial velocity 

CN( ' field were obtained. We analysed the 3D-kinematics by combining the knot radial velocities, 

, derived from IFS data, with the knot proper motions derived from multi-epoch, narrow-band 

images. The intensity maps built from IFS data reproduced well the morphology found in the 
1/-} \ narrow-band images. We checked the results obtained from previous long-slit observations 

. with those derived from IFS spectra extracted with a similar spatial sampling. At the positions 

CS| ' intersected by the slit, the physical conditions and kinematics derived from IFS are compatible 

t-H I with those derived from long-slit data. In contrast, significant discrepancies were found when 

the results from long-slit data were compared with the ones derived from IFS spectra extracted 
• i-H ' at positions shifted a few arcsec from those intersected by the slit. This clearly revealed IFS 

^ \ observations as the best choice to get a reliable picture of the HH emission properties. 

H ■ 

C3 ■ Key words: ISM: jets and outflows - ISM: individual: HH 223, L723, VLA 2 Techniques: 

Imaging spectroscopy: Integral Field Spectroscopy 



1 INTRODUCTION 

The Herbig-Haro object 223 of the Reipurth Catalogue dReipurthl 
1 19941) is a knotty, undulating emission of ~ 30 arcsec length, and 
corresponds to the brightest optical counterpart of a large-scale 
(~ 5.5 arcmin) H2 outflow. The outflow is located i n the dark cloud 
Lynds 723 (L723), at a distance of 300 ± 150 pc dGoldsmith et all 
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Il984|). It is associate d with the low-mass multiple system VLA 2 
dAnglada et al.ll 199lh , being one of the c omponents of this system 
(VLA 2A) the outflow exciting source dCarrasco-Gonzalez et al.1 
2008). For a more detailed description of the region and the nomen- 
clature of the outflow knots, see lLopez et al.l(l2010al) and references 
therein. 

In a previous work dLopez etal.ll2009h . we performed long- 
slit spectroscopy in the spectral range 5800-8300 A along two slit 
positions, which intersect the peaks of the bright HH 223 knots (A 
to E), including emission from the low-brightness nebula surround- 
ing the knots. The spectra of the knots are characteristic of shock- 
excited gas with an intermediate/high degree of excitation. The ra- 
dial velocities derived for the knots are highly blueshifted relative 
to the ambient gas cloud, and supersonic, as expected for a HH op- 
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tical jet. Significant changes in the physical conditions (namely, in 
the excitation, ionization and electron density) were found as we 
move a few arc sec along the slits. In addition, two velocity com- 
ponents seem to be present at several of the positions intersected 
by the slits. It is likely that, at these positions, there is contribution 
to the emission from both the shocked jet gas and the ambient gas, 
dragged and accelerated by the supersonic outflow. Note, however, 
that we sampled partially the spatial emission of the HH 223 knots 
and interknot nebula by means of the long- slit observations. In con- 
trast, Integral Field Spectroscopy (IFS) appears to be the appropri- 
ate tool to get a wider spatial coverage that includes the whole spa- 
tial knot/interknot emission. This technique has proven to be very 
efficient to study the kinematics and physical conditions, mainly 
in the case of Herbig-Haro objects having a p oorly colli mated mor- 
pholog y (see, e. g. HH 262. lL6pez et alJ2008l : HH 1 lO. lLopez et aD 
l2010bh . 

Taking advan tage of the capabilities of INTEGRAL 
dArribas et al.l [l998), we made a single-pointing observation in 
HH 223, covering the emission of its brightest knots (A to D). From 
these data, we obtained the morphology in several emission lines, 
the kinematics and the physical conditions of HH 223. In addition, 
we were able to compare the results derived from IFS data with 
those obtained from narrow-band and long- slit observations using 
similar spatial and spectral resolutions. This experiment gave us 
the opportunity to test whether a partial sampling of the emission 
performed from long- slit data might mask some behaviours of the 
actual spatial distribution of the physical conditions in HHs such 
as HH 223, with a poor collimated and rather chaotic morphology. 

We also present in this paper the 3D kinematics of the bright- 
est HH 223 knots derived by combining the INTEGRAL data (ra- 
dial velocity) with narrow-band Ha images obtained at three dif- 
ferent epochs (proper motions). 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 INTEGRAL FIELD SPECTROSCOPY 

The 4.2 m William Herschel Telescope (WHT) of the Ob- 
servatorio del Roque de los Muchachos (ORM, La Palma, 
Spain) was used in combinatio n with the INTEGRAL fiber op- 
tics system dArribas et al.lll998l) and the WYFFOS spectrograph 
(Bing ham et al.ll 19961) . INTEGRAL links the Nasmyth focus of the 
WHT with the slit of WYFFOS through three optical fiber bundles. 
At the focal plane, the fibers of these three bundles are arranged into 
two groups, one forming a central rectangle, and the other an outer 
ring mapping the sky. The three bundles can be interchanged online 
depending on the scientific program or the prevailing seeing con- 
ditions, as they have different spatial sampling and coverage. More 
information about INTEGRAL bundles and performances can be 
found at the INTEGRAL web pag^E The data discussed in this pa- 
per were obtained with the standard bundle 2 (STD2). STD2 con- 
sists of 219 fibers (189 science and 30 sky fibers), each of 0.9 arc- 
sec in diameter on the sky and covering a field of view (FOV) of 
16 x 12 arcsec 2 . 

The observations were made on 2009 May 25 as a backup 
program of the science verification allocated time for a new equal- 
ized fiber bundle for INTEGRAL. The WYFFOS spectrograph was 
equipped with a 1200 groove mm -1 grating (centered on A 7000 A) 
and a EEV mosaic (two EEV-42-80 thinned and AR coated CCDs) 

1 http :// www. iac . es/proy ecto/integral/ 



of 4308 x 4200 pixels of 13.5 fim in size, giving a linear dispersion 
of about 0.4 A pixel -1 in the spectral range ~ 6180-7700 A. With 
this configuration, and pointing to the brightest knots of HH 223 
(a=\9 h 17 m 57!7, 6=+\9° 11' 52"), five exposures of 1800 s each 
and one exposure of 900 s were taken to get a total effective in- 
tegration time of 2.75 hr. This single pointing covers the brightest 
HH 223 knots (HH 223-A to -D) and the low-brightness nebular 
emission surrounding these knots, as illustrated in Fig.[T] where the 
spectrum extracted from each fiber, in a wavelength range around 
Ha, is overplotted o n the HH 223 Ha narrow-band image. The 
knots are labeled as in Lop ez et al . (200|). 

The data were reduced usins IRAI0 standard routines (see 
e. g. lGarcia-Lorenzo et al.l20 05 for a description of the reduction of 
IFS data from fiber-based instruments). We obtained typical wave- 
length calibration errors of 0.1 A, which give velocity uncertain- 
ties of about ± 4.5 km s _1 for Ha. The differential atmosp heric 
refract ion was estimated according to the model given by lAUenl 
dl973h . being this effect negligible for the HH 223 observations. 
The night was photometric according to the Mercator Telescope 
measurement^! and the seeing was around 1 arcsec. 

2.2 Narrow-band CCD images 

Deep narrow-band CCD images of L723, covering a field of 
~ 5x5 arcmin 2 that includes HH 223 were obtained at three 
different epochs (2004 July 20, 2007 July 14, and 2009 May 18), 
suitable to measure proper motion displacements of the knots. 
The Andalucia Faint Object Spectrograph and Camera (ALFOSC) 
was used on the 2.6 m Nordic Optical Telescope (NOT) at the 
ORM. The image scale was 0.188 arcsec pixel -1 . The filter used 
had a central wavelength of 6564 A and a bandpass of 33 A, 
which mainly included emission from the Ha line, with some 
contamination from the [Nil] A 6548 A emission line. All the 
images were processed as described in Lope z et al 1 (l2006l) for the 
2004 observing run. The IAA Guaranteed Observing Time with 
ALFOSC was used for the observing run of 2009. 



3 IFS RESULTS 

3.1 IFS imaging: spatial distribution of physical conditions 

Emissions from Ha, [N II] and [S II] A 6716, 6731 A were detected 
in all the fibers covering the HH 223 knots and its surrounding 
nebula mapped with INTEGRAL. Other characteristic HH emis- 
sion lines included in the observed wavelength range ([O I] doublet, 
[Fell] A 7155 A and [O II] A 7300 A doublet) were only detected 
in a few fibers, located around the peak positions of the brightest 
knots A and B. 

Figure [2] shows the morphology of HH 223 for the Ha, 
[Nil] A 6584 A and [Sll] A 6716, 6731 A emission lines ob- 
tained from IFS data. To generate these narrow-band images, a 
two-dimensional interpolatio n was applied using the IDA package 
dGarcfa-Lorenzo et al .1120021) . In particular, an ASCII file with the 
actual position of the fiber and the corresponding spectral feature 

2 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 

3 http://www.ast.cam.ac.uk/dwe/SRF/camc-extinction.html 
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Figure 1. (a) Ha narrow-band image of the HH 223 field obtained at the Nordic Optical Telescope (NOT) on 2009 May 18. The white rectangle includes the 
region mapped with the IFS single pointing, (b) Two-dimensional distribution diagram of the observed spectra in the spectral window 6540-6590 A, which 
includes the Ha and [Nil] 6548, 6584 A emission li nes. The spectra have been rescaled to show better the line profile, and are overplotted onto the Ha 
narrow-band. Labels of knots identification (following lLopez et al.ll2006h are included. The offsets are relative to the centre of the field covered by IFS. 
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Figure 2. IFS-derived narrow-band images for several emission lines (Ho;, [Nil] A 6584 A and [S n]A 6716, 6731 A). For each position, the flux of the line 
was obtained from a Gaussian fit to the line profile, after subtracting a continuum, obtained from the adjacent wavelength range free of line emission. The flux 
scale is in counts. 



was transformed to obtain a regularly spaced rectangular grid. In 
this way, we built up narrow-band images with a field of view of 
47x27 pixels (~ 14 arcsec x 8 arcsec) with a spatial scale of 0.3 
arcsec pixel -1 that can be treated with standard astronomical soft- 
ware packages. 

The IFS maps (Fig. [3 are in good agreement w ith previous 
CCD narrow-band Ha and [S II] images of HH 223 dLopez et"aD 
2006). In the four emission lines considered, the IFS maps repro- 
duce the undulating, knotty structure of the emission, surrounded 
by a low-brightness nebula that had already been found in the 
narrow-band, CCD images. Note however that the IFS data allowed 



us to map separately the emission coming from the Ha and [N II] 
lines, as well as the two lines of the [S II] doublet. We can now 
observe the close similarity between the brightness distribution of 
the [S II] and [N II] maps. In contrast, the comparison of Ha and 
[N II] IFS maps shows differences in the morphology and relative 
brightness of the knots. In particular, knot C shows a somewhat dif- 
ferent morphology in [N II] (with two, north- south, condensations 
of similar brightness) and Ha (where the southern condensation is 
fainter than the northern one). Furthermore, knots B and C show a 
higher [N Il]/Ha brightness ratio than knot A. As we will discuss 
later, these signatures are the result of differences in gas excita- 
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Figure 3. [S n]6716/6731 (top), [S n] (671 6+673 1)/Ha (middle) and [N II] 
6584/Hcn (bottom) line-ratio maps obtained from the IFS-narrow-band im- 
ages, blanking the regions with a SNR ^ 3. Contours from the Ha IFS im- 
age are overlapped to help with the identification of the knot condensations. 



tion through the HH 223 emission, which were only outlined from 
long- slit spectroscopy. 

Line-ratio maps of the emission lines widely used for nebular 
diagnosis ([S II] (671 6+673 l)/Ha and [Nil] 6584/Ha for excita- 
tion, and [S II] 6716/6731 for electron density) were obtained from 
the IFS-narrow-band images, and are displayed in Fig. [3] 

For all the mapped region, which included emission from 
knots and from the low-brightness interknot nebula, the values of 
the [Sll]/Ha line ratio range from 0.6 to 1.4, and those of the 
[N ll]/Ha, from 0.2 to 0.5. Accordin g to th e [S ll]/Ha vs. [N ll]/Ha 
diagnostic diagrams (see e. g. Canto 198lh . these values are found 
in HH objects, where the emission is produced by shocks. In fact, 
the values derived for the [S ll]/Ha line-ratio correspond to the 
spectra with an intermediate/high degre e of excitation, fo llowing 
the classification for HHs proposed by iRaga et al.1 dl996h . which 
distinguished the low (with a [S Il]/Ha ^ 1.5) from the intermedi- 
ate/high excitation spectra. The spatial distribution of the electron 
density (n e ) is outlined in the [S II] 6716/6731 line ratio map. Val- 
ues of the sulphur line ratio range from 0.6 to 1 .2, which correspond 
to n e ranging from 0.24 to 4.2 xlO 3 cm -3 (n e derived using the 
TEMDEN task of the IRAF/STSDAS package, and assuming T e =10 4 
K). 

A previous work (Lopez etalJ l2009h partially sampled the 
physical conditions in HH 223 by means of long-slit observations, 
with a slit positioned along HH 223 crossing the knot intensity 
peaks. The long- slit data revealed that knot A had the highest ex- 
citation spectrum. A trend of increasing excitation going towards 
the west was found for the rest of the knots. The n e values derived 
for knots A and B were significant higher (by factors of ~ 5 to 10) 



than for the rest of the knots. Hence, the line-ratio values found for 
the positions observed in the long- slit HH 223 sampling are con- 
sistent with those obtained, at these positions, from the line-ratio 
IFS maps. Note, however, that the IFS data are better suited to map 
the physical conditions in the whole HH 223 region covered by the 
single IFS pointing. As a result, the maps of Fig. [3] show a rather 
more complex spatial distribution for both, excitation and density, 
which was not self-evident from the partial sampling of long-slit 
data. 

The electron density (n e ), derived from the [S II] 6716/6731 
map, is higher in the southern, eastern knots (A, B), with n e ~ 
4200-2700 cm -3 respectively, than in the other two mapped knots 
(C, D), with n e ~ 650-500 cm -3 . In addition, the n e map has 
revealed other density enhancements offset from the knots inten- 
sity peaks (e. g. the one towards the northwest of knot B, within 
the low-brightness nebula connecting knots B and C, with n e ~ 
1400 cm -3 ), as well as a variation of n e across knot B, moving 
from the northwest to the southeast, and ranging from ~ 850 to ^ 
2700 cm -3 . Hence, the [S II] 6716/6731 map likely suggests den- 
sity inhomogeneities within the emitting gas at smaller scales than 
the knot sizes. 

The spatial distribution of the gas excitation is outlined in the 
[S Il]/Ha and [N Il]/Ha maps of Fig.[3] The emission with the high- 
est degree of excitation is found for knot A, and also around the 
peak intensity of knot C, while the emission associated with knot 
B has the lowest excitation degree. However, like for the spatial 
distribution of the electron density, the overall pattern is very com- 
plex, suggesting again smaller- scale inhomogeneities in the excita- 
tion conditions within the mapped region. 

3.2 IFS radial velocity field 

Radial velocity maps were first derived from the line centroids of 
a single-Gaussian fit to the Ha and [S II] emission line profiles. 
However, a careful inspection of the spectra extracted at several 
positions on HH 223 reveals that the emission lines display broad, 
asymmetric and even double-peaked shapes, suggesting that at least 
two kinematic components are contributing to the line profile at 
these positions. Hence, we also derived the radial velocities by ap- 
plying the cross-correlation technique through the XCORR rou- 
tine of the DIPSO package. Ideally, the template must be a high- 
S/N spectrum, with a resolution similar to or better than the tar- 
get spectrum. Unfortunately, we lack such an ideal template, and 
we adopted as template the spectrum extracted at the peak of the 
brightest knot (knot A), which seems to have symmetric line pro- 
files. The velocity maps derived from these two methods are in 
good agreement and show a rather complex pattern in the observed 
lines. Figure |4] (upper panel) shows the maps of the radial veloc- 
ity field for the Ha and [S II] emissions. In order to estimate the 
radial-velocity uncertainties, we compared the values derived from 
the two techniques (single-Gaussian fit to the profiles, and cross- 
correlation) giving a mean difference in radial velocities of 2±6 
km s _1 . The largest differences between the radial velocity maps 
derived from the two techniques are found at positions with clear 
double-peaked line profiles. 

The derived velocitiefl are highly blueshifted, with values 
ranging from -180 to -60 km s _1 . The velocity field shows a rather 

4 Velocities are referred to the local standard of rest (LSR) frame. 
A Vlsr=+10-9 km s _1 for the parent cloud has been taken from 
iTorrelles et aDl986l 
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Figure 4. Lower panels: Integrated intensity maps obtained by integrating the signal within all the Ha (left) and [S n] A 6716 A (right) line profiles (zero- 
order intensity momenta). For each of the panels, the red contours correspond to the emission intensity obtained by integrating the signal for Vlsr ^ -90 
km s _1 (LVC) and the blue contours, to the emission intensity for Vlsr ^ -90 km s _1 . Upper panels: Radial velocity field obtained from cross-correlation 
in the Ha+[N n] range (left), and [S n] 6716, 6731 A ranges (right). The velocity scale is in km s _1 . White contours, as in Fig. [3] 



complex pattern. In spite of that, some trends can be observed in 
the spatial distribution of the velocities: the lowest (absolute) ve- 
locity values, ranging from -80 to -60 km s _1 , come from the 
low-brightness emission that connects knot A with knots C and 
D, and from a shell at the south of the knots A to D intensity 
peaks. The highest (absolute) velocity values, ranging from -180 
to -150 km s _1 , come from a region at the northeast of knot C in- 
tensity peak. Finally, we derived velocities ranging from c± -130 to 
~ -90 km s _1 for the emissions arising around the knot A, C and 
D intensity peaks. 

The line profile shapes of most of the spectra suggested a 
contribution to the emission from two kinematic components. To 
explore this possibility, we built two chanel maps, by consider- 
ing two ranges for the velocity: the higher blueshifted component 
(HVC), for Vlsr^ -90km s _1 , and the lower blueshifted compo- 
nent (LVC), for VisR^ -90km s _1 . Then, we obtained the emis- 
sion maps of the LVC and HVC, for the Ha and [S II] lines, by 
integrating the signal within the corresponding velocity range. Fig- 
ure!?] (lower panels) displays the spatial distribution of the two ve- 
locity components. Contours for the HVC (blue) and LVC (red) 
have been overlapped onto the integrated emission (zero-order in- 
tensity moment). All the knots have emission from both kinematic 
components. However, the spatial distribution of these two velocity 
components does not fully coincide in Ha neither in [S II] lines. 
The most remarkable differences appear for knots B and C, where 
the peak intensity positions of the HVC and LVC appear shifted. 
We measured offset values of 0.9 and 1.5 arcsec in Ha and [S II] 
respectively for knot C, and offsets of 0.7 and 1.5 arcsec for knot B. 
These offset values are reliable, according to the spatial resolution 
of the data. 

Clear doubled-peaked emission lines were only found in a set 
of spectra arising from the bright nebula at several positions south 
of the knot A and B intensity peak, and from the region cover- 
ing knot C (see Fig. [5] upper panel). We additionally performed a 
line-profile decomposition to reproduce the Ha+[Nll] 6548,6584 



A, and the [Sll] 6716,6731 A line profiles of the spectra coming 
from these regions. The line-profile decomposition was obtained 
from a model including two Gaussian components. The fitting was 
performed using the DIPSO package of STARLINK, and the re- 
sults are shown in Fig. [5] The two velocity components derived in 
this way appear blueshifted. For the higher blueshifted component 
(HVC), velocities ranging from -180 km s" 1 to -130 km s" 1 are 
obtained. The velocities derived for the lower blueshifted compo- 
nent (LVC) range from -80 km s _1 to -40 km s _1 . 

Unfortunately, there are few spectra where a reliable line- 
profile decomposition can be obtained, not allowing us to infer 
whether the two kinematic components are tracing emission with 
different physical conditions, having thus different origin or not. 
However, a trend is found that consists in the HVC having a higher 
degree of excitation (with [S Il]/Ha ratios ^ 1) and lower electron 
density (n e ranging from ~ 200 to 800 cm -3 ) than the LVC, which 
is less excited ([Sll]/Ha ^ 1) and denser (n e ranging from 1500 
to 6000 cm -3 ). It should be mentioned that a similar behaviour for 
the excitation and electron density have been detected in jets emerg- 
ing fr om Class 0/1 sources and in microjets from TTauri stars (see, 
e. ^.lHammar]|l994lHrrt. Mundt & Sol3ll99llGarcfa-L6pez etail 
2008; lGarcia-L6pez et al. 2010h . From the analysis of the double- 
peaked line profiles, these authors found that the emission from the 
LVC usually presents higher electron density and lower excitation 
than the emission coming from the HVC. In most cases, the LVC 
is spatially more compact, and is detected only up at a few arcsec 
near the exciting jet source, around the jet launching region, be- 
ing related to the launching mechanisms. This scenario does not fit 
our case, because the knots mapped in HH 223 are far away (~ 
0. 1 pc) from the outflow exciting source. Two-peaked line profiles 
have also been observed i n th e brighter knots o f HH 32 (see e. g. 
ISolf. Bohm & Radl 19861 . and lBeck et aD 1 20041) . In this case, the 
profiles have been interpreted as originated from the emission of 
bow- shocks. However, at the light of the 3D kinematics derived 
from our observations, it is not clear that the bow-shock interpreta- 
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Table 1. Proper motions (Vr) 1 and full spatial velocities (Hot) of HH 223 


knots 










Knot 


Vt 


PA 


Vtot 






(kms- 1 ) 


(degrees) 


(kms" 1 ) 


(degrees) 


A 


36.4 ±2.2 


144 ±6 


107.6 ±5.0 


70 ±5 


B 


49.2 ±3.9 


106 ±3 


84.4 ±6.4 


54 ±6 


C 


3.0 ±6.5 


52 ±104 


120.7 ±8.3 


89 ± 19 


D 


15.0 ±2.4 


99 ±6 


120.0 ±8.7 


89 ± 10 


E 


60.3 ±2.7 


87 ±1 






Fl 


40.8 ±3.9 


80 ±11 






F2 


64.3 ±2.5 


91 ±6 






F3 


7.5 ±4.2 


194 ± 57 







1 A distance of 300 pc has been adopted. The uncertainty in the L723 
distance was not taken into account for derive the errors in proper motions 
quoted in the table. 

tion could work so well in HH 223 (see later). In other cases (e. g. 
HH 34 and HH 46-47) the origin of the LVC at larger distances 
from the source is attributed to the gas entrained and accelerated 
by the high- velocity outflow. This might be a plausible explanation 
for the origin of the LVC emission in HH 223. However, a more 
close similarity in the spatial distribution of the two velocity com- 
ponents should be expected in this case. Alternatively, the detection 
of two velocity components that also differ in their excitation has 
been attri buted to episodic velocity variations in the outflow (e. g. 
HH 15l jMovsessian. Magakian & Moiseevll 20121 ). This possibility 
(i. e. episodic mass-ejection events with variable velocity and direc- 
tion, as can be expected when the exciting outflow source belongs 
to a multiple system) might be a more plausible explanation for the 
origin of the LVC and HVC emissions in HH 223. 



3.3 Long-slit vs. IFS 

In spite of the comparisons not being fully straightforward, one of 
the goals of these observations was to check whether the results de- 
rived from IFS data are fully consistent with the ones obtained from 
the long-slit sampling. Apparently some discrepancies in the spa- 
tial behaviour of the velocity appeared when IFS results (this work) 
are compared with those obtained from long-slit observations of 
lLopez et aD d2009h . For example, as can be seen from Fig. |4] low 
negative velocity values are found associated with knot B emis- 
sion (Vlsr— -70 km s _1 around the Ha knot intensity peak), and 
high negative velocity values appear associated with knot C (with 
Vlsr— -150 km s _1 around the Ha knot intensity peak), reach- 
ing even higher negative values (Vlsr^ -180 km s _1 ) at positions 
north to the knot C intensity peak. In contrast, we found velocities 
of ~ -130 km s _1 and c± -70 km s _1 at positions sampled from 
long-slit spectroscopy through knots B and C, respectively. How- 
ever, we concluded that, in this case, the differences found between 
the results from IFS and long- slit observations mainly arise as a re- 
sult of an incomplete sampling of the highly complex kinematics 
ofHH223. 

To make more meaningful the comparisons, we proceeded 
to simulate long- slit spectra from the IFS data with a close simi- 
lar sampling as in long- slit observations. First, we built a grid of 
0.3 arcsec spacing; then, we interpolated the signal to extract one- 
dimensional spectra at positions intersected by two simulated slits 
of 1 arcsec width. The simulated slits were positioned crossing the 
peak intensities of knots A and C, at position angles (PA) 73° and 
79° respectively. We derived the radial velocities (from the cen- 
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Figure 5. Top: Location of the spectra where double-peaked emission lines 
are detected, overplotted onto the Ha narrow-band image of HH 223 (same 
as Fig[T}. Bottom: Observed spectra (black, continuous) and fits (blue, 
dashed) obtained from the model of two Gaussian components applied to 
the Ha+[Nn] 6548,6584 A (left) and [Sn] 6716,6731 A (right). The two 
individual components are shown below in each panel. The same colour 
(red and blue) indicates the same component for each line. Residual ob- 
tained subtracting the fit from the observed spectra are plotted in green. 
Labels in the right corner of each plot indicate the corresponding spectra 
plotted in the top panel of this figure. 



troids) and line ratios (from the fluxes) of the HH 223 emission 
within the simulated slits from Gaussian fits to the lines of these 
extracted spectra. Figure |6]shows the radial velocities derived from 
the centroids of the Ha and [Sll] 6716 A emission lines (blue). 
The corresponding radial velocities derived from long-slit spectra 
are represented by the black dots. Figure |7]shows the [N ll]/Ha and 
[S Il]/Ha line ratios (both are tracers of the excitation conditions) 
and the electron density, n e , from the [S II] 6716/6731 line ratio, 
obtained from the spectra along the simulated slit (blue) and the 
corresponding values from long-slit observations (black). As can be 
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Offset(arcsec) 

Figure 6. Comparison of radial velocities derived from long-slit spectra 
(black) and from IFS spectra, obtained through simulated slits of 1 arcsec 
width positioned at P.A. 79°: crossing the knot C intensity peak (blue); 
the same, but 2 arcsec north of it (red). Offsets are relative to the knot C 
intensity peak. 




10 

X offset (arcsec) 

Figure 8. CCD Ha image of HH 223 from the last epoch (2009). Arrows 
indicate the proper motion of the knot. The ellipse at the end of the arrow 
indicates the uncertainty in the components of the velocity vector. Boxes 
include the regions used for calculation of the knot proper motion. The same 
boxes have been taken to evaluate the representative radial velocity of the 
knot used to derived the full spatial velocity (see text). 
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Figure 7. Comparison of line ratios obtained from long-slit spectra (black) 
and from IFS spectra obtained through simulated slits of 1 arcsec width and 
positioned at P.A. 79° crossing the knot C intensity peak (blue). 



seen from these figures, both the kinematics and the physical con- 
ditions of the emission intersected by the simulated slit are com- 
patible with those found from the long-slit sampling. In contrast, 
some discrepancies appear when velocities derived from long- slit 
are compared with those obtained from IFS spectra through the 
simulated slit, but displaced, in parallel, ~ 2 arcsec north to the 
one crossing the knots peaks (red dots in Fig. [6}. This test reveals 
the importance of obtaining 2D maps covering the full spatial ex- 
tension of the emission, instead of doing ID long-slit sampling, 
mostly in the cases where the kinematics and physical conditions 
of the emission show such as complex pattern, as in HH 223. 



tions. With this aim, we combined the proper motions of HH 223, 
derived for the first time in this work, and radial velocities, derived 
from the current IFS data. 



4.1 Proper motion determination 

The tangential velocity, Vt, of the knots was derived from their 
proper motions, using three narrow-band Ha images of the L723 
field, obtained at the NOT telescope, with a time baseline of 5 years 
(see §2.2). These images were converted onto a common refer- 
ence system using the position of ten common field stars, well dis- 
tributed around HH 223, to register the images. The GEOMAP and 
GEOTRAN task of IRAF were applied to perform a linear transfor- 
mation with six free parameters that take into account relative trans- 
lation, rotation and magnification among the frames. The procedure 
to calculate the prop er motions from mu ltie poch images has been 
described in detail by Lopez et al 1 d2QQ5h and lAnglada et al ] j2007h 
for the cases of HH 110 and HH 30, respectively. In short, it con- 
sists of defining boxes in each frame that include the emission from 
the individual knot; then, to compute the two-dimensional cross- 
correlation function of the emission for the two pairs of frames (the 
2004 first epoch image and each of the other two epoch images), 
and to determine the displacements in x and y coordinates through 
a parabolic fit to the peak of the cross-correlation function; finally, 
to perform a linear regression fit to the displacements of the x and 
y coordinates as a function of the epoch offset, from the first epoch 
(2004). We derived the proper motion velocity (Vt) and the posi- 
tion angle (PA) for each knot from the slope of the fitted straight 
line (and its uncertainty) . The results are displayed in Fig. [8] and 
their values, in columns 2 and 3 of Table Q] The errors include the 
uncertaint y in the position o f the correlation peak (derived as ex- 
plained in lLopez et "all l2005h . and the error in the alignment of the 
frames. 



4 3D KINEMATICS 

The HH 223 emission shows an appreciable undulating morphol- 
ogy. To better characterize the velocity field of the emission, we 
searched for the full spatial velocity and the inclination angle in the 
plane of the sky for the HH 223 knots covered by the IFS observa- 



4.2 Full spatial velocities 

The IFS data have revealed a complex kinematic structure for the 
HH 223 emission in such a way that radial velocity values of the 
knots show variations depending on the parcel of the knot emission 
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Table 2. V LS r (km s _1 ) line centroids 1 of the HH 223 knots from IFS data. 



Knot [Nil] 


Ha 


[Nil] 


[Sii] 


[S II] Average rms 


6548 A 




6583 A 


6716 A 


6731 A 



A 


-88.3 


-101.3 


-86.9 


-86.1 


-85.4 


-89.6 


6.6 


B 


-78.7 


-68.6 


-77.3 


-67.8 


-67.2 


-71.9 


5.6 


C 


-94.2 


-120.7 


-91.9 


-103.0 


-102.3 


-102.5 


11.3 


D 


-116.7 


-124.3 


-115.2 


-120.7 


-120.0 


-119.4 


3.6 



1 Typical value of the error is 5 km s 1 



sampled within the slit (see e. g. Fig. [6}. Note that the emission of 
each knot in HH 223 extends over several arcsec, both in the x and 
y coordinates (see e. g. Fig. [TJ, being thus improperly sampled by 
a single long-slit ponting. Since IFS allows us to achieve a more 
complete sampling of the kinematics, we can assign to each knot 
a representative radial velocity coming from the emission enclosed 
in the same box that we defined to measure the knot proper motion. 
The radial velocity derived in this way will be used to calculate the 
full spatial velocity of the knot. We then defined the aperture en- 
closing each of the boxes A to D (see Fig. [8]), and compute the line 
profiles by integrating the signal over the fibers inside the aperture. 
We then determined the radial velocity from a Gaussian fit to the 
emission lines (see Table 0. 

The full spatial velocity (Vtot) and the angle between the 
knot motion and the plane of the sky (with positive values of <j) to- 
wards the observer) are then derived from the Ha radial velocity 
(obtained from IFS data) and the proper motions (obtained from 
Ha CCD narrow-band images). The results are given in columns 
3 and 4 of Table Q] Unfortunately, we were only able to obtain Vt 
and cj) for the four HH 223 knots covered by the single IFS point- 
ing. We found that these knots present high full spatial velocities 
(^ 100 km s" 1 ) and move in a direction close perpendicular to 
the plane of the sky. However, the data suggest that there is some 
change in the jet inclination, which increases going from the south- 
ern knot B (also with the lowest full velocity) to the northern knots 
C and D. 

Apparently, the 3D kinematics of HH 223 (i. e. the radial ve- 
locity structure and th e proper motions) shows close similarities 
with HH 32 (see, e. g. lSolf. Bohm & Ragd[l98fl and lCuriel etail 
Il997h : both HHs are oriented close to the line of sight and present 
two-peaked line profiles (note however that in the case of HH 223, 
the two-peak ed line profiles ar e not detected at the brightest knot, 
HH 223-A). iRaga et alJ <2004h modeled the HH 32 emission as 
a broken-up, large scale bow- shock with several condensations. 
The bow- shock interpretation succeeded in deriving the bow -shock 
properties, like the bow- shock velocity (Vbs) and the orientation 
angle of the outflow in the plane of the sky (0), from the com- 
ponents of the proper motion velocity paralle l and perpendicular 
to the projected outflow axis dRaga et al.ll 19970 . However, in spite 
of the similarities between the HH 32 and 223 structures, the bow- 
shock interpretation seems not work so well for HH 223. Assuming 
a projected outflow axis at RA. of 110-115° (the RA. of the east- 
west CO bipola r outflow, see e. g. lMoriartv-Schieven & Snell l989: 
iLee et al.l F2002) and from the measured proper motions, the bow- 
shock model would give a low shock velocity (Vt> s ^ 25 km s _1 ) 
and an orientation outflow angle close to the plane of the sky {<j> ^ 
30°), which is far from the direction derived using the observed ra- 
dial and tangential velocity components presented here, even taking 
into account the uncertainty of the proper motions due to the uncer- 
tainty in the distance to L723. 



5 SUMMARY AND CONCLUSIONS 

We carried out a single-pointing IFS observation of HH 223, the 
optical counterpart of a larger scale (~ 5.5 arcmin, ^ 0.5 pc at a 
distance of 300 pc) H2 outflow, driven by the protostellar source 
VLA 2 A, in L723. Because of its peculiar and poorly collimated 
morphology, this target was selected for the backup list of a science 
verification program of a new INTEGRAL equalized fiber bundle. 
IFS observations have revealed to be the best choice to achieve 
a reliable mapping of the whole spatial emission of such chaotic 
outflows. The results derived from IFS data were compared with 
previous long- slit and CCD narrow-band images looking for con- 
sistency. Unfortunately, because of the observing time available for 
our target in the backup program, it was not possible to acquire an 
additional pointing to cover the fainter HH 223 knots E to F. New 
results derived from our analysis of the IFS data are summarized 
below: 

• The morphology of the HH 223 emission shown by the line 
intensity maps built from the IFS data is in good agreement with 
the morphology found in the narrow-band images. One advantage 
to get line intensity maps from IFS data is the capability of mapping 
separately the emission coming from lines that usually are included 
into the same narrow-band filter (e. g. Ha and [Nil] 6584 A). We 
obtained in this work a map for each of the emission lines that 
revealed the close similarity between the [S 11] and [N 11] emissions 
and some differences between the spatial brightness distribution of 
the neutral (Ha) and ionized ([Nil] and [S II]) gas. 

• We traced the spatial distribution of the excitation and the 
electron density from line-ratio maps. These maps revealed the 
complexity of the spatial structure of the physical conditions. It 
makes once again IFS mapping highly appropriate to properly re- 
veal physical conditions of the emission, instead of derive them 
from a partial sampling obtained from long- slit data. As a gen- 
eral trend, the [S ll]/Ha and [N ll]/Ha ratios are consistent with 
the emission arising from shocks with an intermediate/high degree 
of excitation. A relevant result concerning the electron density is 
the inhomogeneities found through the emission (both, through the 
knots and through the low-brightness nebula), which suggests a 
clumpy structure that cannot be well resolved with the current IFS 
spatial resolution. 

• The radial velocity field has been derived, first from the line 
centroid of a single-Gaussian fitting to the lines profile and, in ad- 
dition, from the cross-correlation technique. Results found from 
both procedures are fully consistent. The velocities are highly 
blueshifted, ranging from -180 to -60 km s _1 . They present a com- 
plex pattern in the spatial distribution, with some trend that consists 
of the more negative (blueshifted) values increasing moving from 
southeast to northwest. 

• The spectra extracted at several positions show clear double- 
peaked line profiles. We were able to reproduce the line profiles 
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for these spectra using a model that included two blueshifted kine- 
matic components. In addition, some trend is found that consists of 
the HVC (velocity ranging from -180 to -130 km s _1 ) being more 
excited and rarefied than the LVC (velocity ranging from -80 to 
-40 km s _1 ). The two velocity components, having quite different 
physical properties and spatial distributions, can be indicative of a 
variability in the properties of the outflow ejection (e. g. episodic 
mass ejection events having some different velocity and direction). 
This can be expected when the o utflow source is a binary or a mu l- 
tiple system, as was proposed by Carrasco-Gonzale z et al ]{2QQi). 

• To search for the reliability of two kinematic components 
through all the emission mapped, we performed chanel maps both, 
in Ha and in [S II], for two velocity ranges: a HVC (for Vlsr ^ - 
90 km s" 1 ), and a LVC (for Vlsr > -90 km s" 1 ). We found differ- 
ences between the spatial distributions of the HCV and LVC. The 
most remarkable differences are found in knots C and B, where we 
measured offset values from 0.7 to 1.5 arcsec between the posi- 
tions of their LVC and HVC emission peaks. Taking into account 
the spatial resolution of our data, these offsets are reliable. 

• We checked the results on the kinematics and physical con- 
ditions derived from long-slit data with those derived from IFS. 
With this aim, we simulated a long- slit and extracted spectra with 
a similar sampling along this slit, positioned as in long- slit data ac- 
quisition. This test lead us to conclude that the physical conditions 
and kinematics derived from long-slit spectra are in good agree- 
ment with those derived from IFS spectra that were extracted at the 
same positions. In contrast, we found that appreciable differences 
could arise when long- slit data are compared with IFS spectra that 
were extracted along slits displaced few arcsec from the nominal 
long- slit position. This test reinforces again the IFS mode to get a 
more accurate picture of such as complex outflows. 

• Comparing the spatial distribution of the gas excitation (FigO 
and kinematics (FigO, we can observe a similar structure of the 
maps, which illustrates a relationship between the degree of gas 
excitation and the velocity. In particular, the regions having a higher 
degree of gas excitation (around knots C and A) coincide with the 
ones were the velocity values are more negative, while emission of 
lower excitation (located north of knot B and south of knots C and 
D) also displays less negative velocity values. We cannot appreciate 
a similar relationship when the maps of the electron density and 
velocity are compared, probably because the spatial distribution of 
the electron density has a pattern more complex than the excitation 

• We calculated the full spatial velocity and the inclination angle 
of its motion in the plane of the sky for the knots mapped with IFS. 
From proper motion measurements, using multiepoch narrow- band 
images, and radial velocities, derived from IFS, we found that the 
knots have full spatial velocities higher than 100 km s _1 and move 
close to the line of sight. Furthermore, a change in the jet inclination 
is also outlined. 
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